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FOREWORD 


This  study  was  conducted  under  SSD  7993  -  17533,  Search  and  Re¬ 
covery  Concepts,  dated  17  August  I960.  This  task  called  for  a  study  of 
the  problem  of  search  and  recovery  of  both  space  vehicles  and  downed 
aircraft  on  a  global  basis.  Another  task,  SR  7993  -  49759,  Retrieval 
Concepts,  covers  similar  requirements,  and  reports  on  various  aspects 
of  this  task  have  been  prepared  by  several  contractors. 

The  author  wishes  to  acknowledge  the  assistance  of  Mr.  Harry 
Gaines  who  programmed  this  model  for  the  IBM  7090  and  was  of  con¬ 
stant  aid  in  formulating  the  problem.  Also,  Mr.  Gaines  wrote  much  of 
those  portions  of  the  report  dealing  with  the  actual  computer  program. 
The  aid  and  advice  of  Dr.  John  J.  Murphy,  who  was  the  original  project 
officer,  is  also  gratefully  acknowledged. 


APGC-TDR-62-44 


Catalog  catda  with  abatraata  may  be  found  at  the  back  of  thie  document. 


ABSTRACT 


The  subject  of  computer  simulation  was  chosen  for  this  feasibility 
study  as  this  aspect  does  not  appear  to  have  been  well  covered  during 
previous  studies  of  the  Global  Search  and  Recovery  problem. 

A  mathematical  model  was  constructed  of  a  simulated  search  and 
recovery  mission,  and  a  computer  program  was  constructed  which  em¬ 
bodied  this  model.  The  objective  of  this  computer  program  was  the 
determination  of  the  probability  of  successfully  retrieving  a  downed  ob¬ 
ject  when  given  information  as  to  the  probable  location,  the  nature  and 
deployment  of  search  and  recovery  facilities,  and  a  geographic  and 
political  description  of  the  earth.  Considerable  attention  was  given  to 
developing  a  convenient  means  of  entering  and  varying  the  parameters 
of  the  model  in  a  systematic  manner. 

The  results  of  this  study  indicate  that  it  would  be  feasible  to  con¬ 
struct  a  computer  model  of  at  least  certain  aspects  of  the  search  and 
recovery  problem  which  wo\Ud  be  suitable  for  studies  of  the  effects  of 
several  parameters.  However,  the  model  described  is  believed  to  be 
too  limited  and  inflexible  in  form  (though  not  in  parameter  values)  to 
permit  generalisation  to  a  larger  variety  of  search  and  recovery  prob¬ 
lems.  Future  work  along  this  line  should  probably  employ  a  Monte 
Carlo  type  approach. 
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SECTION  1  -  INTRODUCTION 


The  Investigation  of  one  aspect  of  the  search  and  recovery  task, 
that  of  a  computer  model,  was  not  attempted  during  previous  studies. 
Since  this  aspect  was  suggested  as  a  desired  goal  by  Air  Rescue  Serv> 
ice,  such  a  computer  model  was  made  the  subject  of  the  study  covered 
in  this  report.  In  devising  the  model  described,  full  use  was  made  of 
reports  in  the  retrieval  concepts  area  (SR  49759)  made  to  the  Aeronau¬ 
tical  Systems  Division  by  contractors.  Without  such  information,  the 
study  covered  in  this  report  would  have  had  little  relationship  to  the 
actual  problem. 

Since  no  previous  models  applicable  to  this  problem  could  be  found, 
the  initial  phases  of  this  study  were  devoted  to  determining  the  various 
aspects  of  the  problem  and  their  application  to  a  high  speed  computer 
program.  Considerable  thought  was  devoted  to  determining  the  vari¬ 
ables  involved  and  to  how  they  might  be  numerically  specified.  In  par¬ 
ticular,  the  problem  of  specifying  the  pertinent  variables  for  the  surface 
of  the  earth,  on  a  global  scale,  was  most  difficult.  Also  various  basic 
assumptions  had  to  be  made,  for  example,  the  viability  of  a  person  or 
object  was  assumed,  consistent  with  previous  limited  information,  to  be 
a  negative  exponential,  with  one  exponent  up  to  time  of  finding  and  an¬ 
other  exponent  from  time  of  finding  to  time  of  retrieval.  Such  assump¬ 
tions  cannot  be  considered  completely  accurate,  but  it  is  believed  that 
they  are  reasonable  and  not  grossly  incorrect. 

The  specific  search  and  recovery  tactics  to  be  employed  also  cannot 
be  positively  determined;  however,  these  need  be  only  approximations 
having  the  essential  characteristics  of  an  actual  search  and  recovery. 
The  search  tactics  are  considered  to  consist  of  a  dash  by  the  search  ve¬ 
hicle  directly  to  the  most  probable  location  of  the  object,  and  then  a 
spiral  search  about  this  point.  Similarly,  the  retrieval  tactics  are  con¬ 
sidered  to  consist  of  a  dash  by  the  retrieval  vehicle  (possibly  a  ship) 
toward  the  most  probable  position,  then  a  change  of  course,  at  the  time 
the  object  is  sighted,  toward  the  point  at  which  retrieval  can  be  accom¬ 
plished. 

The  computer  model  which  will  be  described  represents  an  attempt 
to  introduce  the  variables  in  a  quantitative  manner  so  that  a  quantitative 
value  of  the  probability  of  retrieval  for  any  given  situation  can  be  found. 
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and  more  important,  the  effects  of  variations  in  the  situation  can  be 
conveniently  studied. 

In  order  to  study  the  problems  involved  in  mechanizing  the  model, 
the  programming  work  was  actually  carried  out.  In  addition,  consider¬ 
able  effort  was  spent  in  investigating  the  possibilities  and  difficulties 
of  varying  almost  any  of  the  parameters  in  a  factorial  manner. 


SECTION  2  -  DESCRIPTION  OF  PARAMETERS  AND  VARIATIONS  TO 
PARAMETERS 


FACTORS  WHICH  INFLUENCE  THE  PROBABILITY  OF  SUCCESSFUL 
RESCUE 

The  study  of  previous  efforts  in  this  field  revealed  that  there  are  a 
number  of  factors  which  are  variable  and  which  greatly  influence  the 
probability  of  a  successful  retrieval.  Some  of  the  most  important  are 
described  in  the  following  paragraphs.  In  the  remainder  of  this  report, 
we  will  refer  to  specific  parameters.  It  should  be  kept  in  mind,  however, 
that  these  parameters  are  formalized  embodiments  of  the  variables  dis¬ 
cussed  below,  and  that  one  of  these  parameters  may  represent  the  com¬ 
bined  effects  of  several  physical  agencies,  devices,  or  variables  of  the 
real  problem.  Thus,  if  a  given  parameter  shows  up  as  crucial,  the 
means  of  obtaining  a  better  value  for  this  parameter  may  require  further 
study.  A  model  such  as  that  described  here  can  only  point  to  the  crucial 
parameters  for  a  given  problem. 

SI  ZE  OF  AREA  TO  BE  SEARCHED.  The  size  of  the  area  to  be 
searched  may  vary  from  a  few  square  miles  to  tens  of  thousands  of 
square  miles.  The  importance  of  this  factor  is  obvious,  but  quantitive 
studies  are  considered  to  be  helpful. 

WIDTH  OF  PATH  WHICH  MAY  BE  SEARCHED  BY  SEARCH  VEHICLE. 
The  width  of  the  path  which  may  be  searched  by  the  search  vehicle  may 
vary  from  a  few  feet  to  a  few  miles  for  a  visual  search,  as  to  a  large  ex¬ 
tent,  visual  search  is  affected  by  weather  conditions.  If  electronic  de¬ 
vices  are  used  with  a  cooperating  object,  the  width  of  the  path  may  be  as 
great  as  several  hundred  miles. 

RANGE  AND  ENDURANCE  OF  SEARCH  AND  RECOVERY  VEHICLES. 
Range  and  endurance  are  controlling  factors  for  long  range  searches,  and 
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design  specifications  for  vehicles  depend  on  the  values  required.  For 
example,  helicopters  usually  cannot  be  used  for  long  range  searches, 
while  ships  perform  well  with  regards  to  range  and  endurance. 

TYPE  OF  RECOVERY  VEHICLE.  The  possibility  exists  that  re¬ 
covery  can  be  accomplished  by  the  search  vehicle,  for  example,  both 
search  and  recovery  may  be  possible  by  an  aircraft  which  can  hover. 

In  other  cases,  separate  vehicles  may  be  required  for  search  and  re¬ 
covery.  This  problem  has  led  to  many  differing  proposals,  such  as 
the  "Mother -Daughter"  system,  large  VERTOL  search  aircraft  and  pick¬ 
up  by  "Long  Line"  techniques. 

POUTICAL  SITUATION  IN  SEARCHED  AREA.  The  political  situ¬ 
ation  in  the  area  of  probable  location,  as  it  affects  the  search  and  the 
possibility  of  retrieval,  is  a  variable  which  realism  requires  to  be  con¬ 
sidered. 

PROBABILITY  OF  SAFE  LANDING.  The  probability  of  the  downed 
object  landing  safely  in  the  area  enters  into  a  decision  for  the  type  of 
retrieval. 

EFFECTS  OF  CLIMATIC  OR  GEOGRAPHICAL  CONDITIONS  ON 
VIABILITY.  The  severity  of  climatic  or  geographical  conditions  in  the 
region,  as  it  affects  the  chance  of  living  until  recovery,  is  an  important 
factor  in  plauining  the  rescue.  Clearly,  a  person  overboard  in  the  North 
Atlantic  must  be  rescued  sooner  than  one  in  more  moderate  tempera¬ 
ture  regions. 

TIME  REQUIRED  TO  LOCATE  AND  RETRIEVE.  The  time  required 
to  locate  and  retrieve  a  downed  object  or  person  has  a  major  influence 
on  the  success  of  a  retrieval,  especially  when  climatis,  or  geographical 
conditions  are  severe.  Such  factors  as  speed  and  scramble  time  of 
search  and  recovery  vehicles  play  a  part  in  determining  the  time  of 
finding  and  retrieval.  The  location  of  bases  and  the  characteristics  of 
the  search  and  recovery  vehicles  stationed  at  each  base  must  also  be 
considered. 


PARAMETERS  AND  TABLES 

The  number  of  parameter  values  necessary  to  describe  the  problem 
configuration  by  use  of  this  model  is  very  large,  approximately  25,000. 

It  was  clear,  therefore,  that  some  convenient  method  was  necessary  to 
introduce  these  values  into  the  computer  and  vary  them  during  an  experi¬ 
ment.  The  parameters  fall  rather  naturally  into  a  small  number  of 
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groups.  These  groups  were  arranged  in  tables,  and  the  tables  were 
arranged  in  a  systematic  way  and  given  names  descriptive  of  the  para¬ 
meters  they  contain. 

The  major  groups  of  parameters  are  those  that  describe: 

1.  The  topographical  and  political  features  of  the  surface  of 
the  earth. 

2.  The  nature  and  deployment  of  search  and  recovery  facili¬ 
ties. 

3.  The  probable  location  of  the  object  of  search  and  recovery. 

4.  The  method  of  search  and  recovery  to  be  used. 

TOPOGRAPHICAL  AND  POLITICAL  FEATURES.  The  group  con¬ 
taining  the  topographical  and  political  features  consists  of  three  tables, 
GLOBE,  PROB,  and  VIAB.  These  tables  are  referred  to  as  permanent 
tables  since  they  contain  information  which  is  less  variable  than  that  con¬ 
tained  in  the  other  tables. 

GLOBE  Table.  The  GLOBE  table  contains  the  description  of  the 
surface  of  the  earth.  The  earth's  surface  was  divided  into  sectors 
bounded  on  the  east  and  west  by  parallels  of  longitude  which  were  4  de¬ 
grees  apart  and  on  the  north  and  south  by  circles  of  latitude  chosen  so 
that  all  the  sectors  had  an  equal  area.  Forty-four  circles  of  latitude 
were  used  together  with  the  poles  so  that  there  were  4050  (45  x  90)  sec¬ 
tors.  These  sectors  are  indexed  with  two  indices,  i  and  j.  The  ij^^  sec¬ 
tor  is  the  i^^  sector  south  from  the  North  Pole  and  the  j^^  sector  east 
from  the  Greenwich  meridian. 

Six  quantities  were  used  to  describe  each  sector: 

1.  Topography. 

2.  Country. 

3.  Probability  of  safe  landing. 

4.  Political  probability  of  retrieval. 

5.  Viability  up  to  time  of  finding. 

6.  Viability  from  finding  to  retrieval. 
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Item  1  was  Included  to  enable  statistics  on  these  features  of  the 
terrain  to  be  obtained. 

Item  2  was  intended  to  provide  a  means  of  changing  item  4  for  all 
sectors  belonging  to  a  given  country.  This  item  is  not  in  the  program  as 
it  now  stands. 

Items  3  through  6  were  used  in  the  calculation  of  the  probabilities. 

Thus,  when  using  this  model,  24,  300  (6  x  40  50)  parameter  values 
are  needed  to  describe  the  surface  of  the  earth.  Since  only  32,  768  words 
of  storage  are  available  in  memory,  this  information  is  packed,  that  is, 
all  six  items  pertaining  to  a  sector  are  stored  in  one  word  and  there  is 
one  word  for  each  sector.  Because  this  packing  limits  the  space  avail- 
able  for  these  items,  most  are  not  stored  directly,  but  rather  an  index 
referring  to  them  is  stored.  Thus,  if  item  5  is  14,  this  means  that  the 
first  viability  parameter  for  that  sector  is  the  14^^  value  in  the  list  of 
viability  values. 

PROB  Table.  PROB  is  a  table  of  32  selected  values  to  be  used  as 
probabilities.  Items  3  and  4  of  the  GLOBE  word  are  indices  referring 
to  PROB. 

VlAB  Table.  VIAB  is  a  table  of  128  selected  values  to  be  used  as 
viability  parameters.  Items  5  and  6  of  the  GLOBE  table  are  indices  re¬ 
ferring  to  VIAB, 

SEARCH  AND  RECOVERY  FACILITIES.  The  parameters  of  the 
second  group,  those  describing  the  search  and  recovery  facilities,  are 
contained  in  the  SBASE  and  RBASE  tables. 

The  SBASE  table  contains  the  description  of  up  to  50  search  bases. 
For  each  search  base,  the  following  information  is  specified: 

1.  Number  of  vehicles  available, 

2.  Velocity  of  vehicles. 

3.  Latitude  of  base, 

4.  Longitude  of  base. 

5.  Scramble  time  of  vehicles, 

6.  Endurance  time  of  vehicles. 
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7.  Probability  of  detection  of  the  object. 

8.  Sweep  width  of  vehicle. 

These  parameters  are  double  indexed,  the  ij^^  entry  in  SBASE  table 
being  the  item  referring  to  the  i*^  base. 

The  RBASE  table  contains  the  description  of  the  retrieval  bases  and 
is  identical  in  format  to  SBASE,  with  the  exception  that  items  7  and  8  do 
not  appear  since  they  are  inapplicable  to  a  retrieval  vehicle. 

PROBABLE  LOCATION  OF  OBJECT,  The  parameters  of  the  third 
group,  those  describing  the  probable  location  of  the  object,  are  stored 
in  the  OBJ  table.  This  table  has  only  three  entries: 

1.  The  latitude  of  the  mean  point  of  impact  (MPI). 

2.  The  longitude  of  the  MPI. 

3.  The  value  of  the  circular  probable  error  (CEP), 

METHOD  OF  SEARCH  AND  RECOVERY.  The  parameters  of  the 
last  group,  those  used  to  designate  the  method  of  search  and  recovery 
to  be  used,  consist  simply  of  two  numbers.  These  two  numbers  are 
kept  in  two  storage  cells  called  respectively  SDOCT,  for  search  doctrine, 
and  RDOCT,  for  retrieval  doctrine.  For  uniformity  of  nomenclature, 
these  are  also  called  tables. 

The  only  flexibility  in  the  search  doctrine  is  the  manner  in  which 
the  search  is  selected.  This  can  be  done  in  one  of  four  ways: 

1,  Choose  the  base  nearest  the  MPI. 

2.  Choose  the  base  whose  vehicle  can  first  arrive  at  the  MPI. 

3.  Choose  the  base  whose  vehicle  can  search  longest  after 
arrival  at  the  MPI. 

4,  Choose  the  base  whose  vehicle  can  search  the  largest  area 
after  arrival  at  the  MPI, 

The  code  number  stored  in  SDOCT  will  determine  which  of  the  above 
methods  will  be  used. 

The  only  flexibility  in  the  doctrine  of  retrieval  is  the  choice  of  whether: 
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1,  The  search  vehicle  also  retrieves;  therefore,  no  retrieval 
bases  are  involved, 

2.  The  search  vehicle  does  not  retrieve;  therefore,  a  re¬ 
trieval  base  must  be  chosen  (always  the  closest). 

The  code  number  stored  in  RDOCT  determines  which  retrieval  doctrine 
will  be  used. 

Note:  When  the  actual  computer  program  is  described,  it  will  be 
found  that  for  technical  reasons  certain  of  the  tables  are  duplicated. 
This  need  not  be  considered  at  this  time  as  this  duplication  was  neces¬ 
sary  only  to  permit  factorial  experiments  to  be  run. 


CALCULATION  OF  A  CASE 

A  single  case,  or  simulated  mission,  consists  of  the  calculation  of 
the  probability  of  successfully  retrieving  a  downed  object,  given  infor¬ 
mation  as  to  the  probable  location,  the  nature  and  deployment  of  search 
and  retrieval  facilities,  and  the  geographic  and  political  conditions  of  the 
area. 


The  method  of  performing  this  calculation  is  to  generate  a  series  of 
representative  locations  (points)  where  the  object  might  be  and  to  calcu¬ 
late  the  probability  of  its  being  successfully  retrieved  from  each  of  these 
points.  A  weighted  average  of  the  probabilities  is  then  calculated,  the 
weights  being  taken  as  the  probability  that  the  object  is  actually  in  the 
area  represented  by  the  points.  The  probability  density  function  of  the 
object  sought  is  taken  to  be  circular  binormal,  the  MPl  and  CEP  being 
furnished.  The  representative  points  are  generated  so  as  to  be  equally 
spaced  in  angle  on  concentric  circles  centered  about  the  MPL  Eight 
points  are  generated  for  each  circle. 


INITIAL  INTRODUCTION  OF  PARAMETER  VALUES 

The  parameters  are  introduced  by  punched  cards  called  SETUP 
cards.  These  cards  are  identified  by  the  word  SETUP  in  columns  1 
through  5  and  contain  in  general: 

1.  The  name  of  the  table  to  be  set  up  beginning  in  column  13. 

2.  Indices  of  flags  denoting  which  part  of  the  table  is  to  be  set 
up. 
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3.  The  value  or  values  to  be  inserted. 

For  details  of  the  format  of  these  cards  see  Appendix  II. 

The  parameter  values  so  introduced  are  the  original  values  for  each 
experiment.  After  being  introduced,  they  are  subject  to  variation. 


VARIATION  OF  PARAMETERS  FOR  A  GROUP  OF  CASES  (OR  EXPERI¬ 
MENT) 

It  is  believed  that  a  "factorial  experiment"  is  a  desirable  and  con¬ 
venient  way  of  systematically  investigating  the  separate  as  well  as  com¬ 
bined  effect  of  a  whole  pattern  of  changes  in  pertinent  parameters.  Such 
factorial  variation  of  parameters  is  specified  by  use  of  FACTOR  cards. 
The  use  of  this  feature  of  the  program  is  not  mandatory,  and  in  the  event 
only  one  case  is  desired  and  the  original  values  are  to  be  as  set  up,  the 
only  card  required  is  a  FACTOR  RESET  card.  Otherwise  a  number  of 
FACTOR  cards  ntuiy  be  employed. 

See  Appendix  I  for  a  discussion  of  a  factorial  experiment.  The  for¬ 
mat  of  FACTOR  cards  is  described  in  Appendix  II. 


CALLING  FOR  CALCULATION  OF  AN  EXPERIMENT 

A  card,  punched  GOGO  in  columns  1  through  4,  is  used  as  a  signal 
that  the  setup  of  an  experiment  is  complete  and  the  calculation  is  to  begin. 

To  summarize,  an  experiment  is  performed  by: 

1.  Entering  a  basic  set  of  parameters  by  means  of  SET  and/or 
LOAD  cards. 

2.  Entering  information  as  to  how  this  basic  set  is  to  be  varied 
during  the  experiment  by  means  of  FACTOR  cards  (this  includes  no  vari¬ 
ation,  which  is  specified  by  a  FACTOR  RESET  card). 

3.  Calling  for  execution  by  means  of  a  GOGO  card. 

This  results  in  calculation  of  all  cases  specified  by  the  factorial  experi¬ 
ment  setup.  Further  experiments  may  now,  if  desired,  be  set  up  and 
run  during  the  same  machine  run. 
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OUTPUT 

For  a  sample  output  see  Appendix  II. 


SECTION  3  -  MATHEMATICS  OF  THE  MODEL 

f 

OVERALL  EQUATION  FOR  PROBABILITY  OF  RETRIEVAL 
Let 


P  s  overall  probability  of  successful  retrieval  for  a  mission. 


Then 


P  =  2,  2  AP(N,  M) 

N  M 

where 

AP(N,  M)  =  incremental  probability  of  successful  retrieval  from  the 

area  (segment  of  a  ring)  having  its  center  at  the  M^^  angu¬ 
lar  position  (©  s  ©i^)  on  the  Nt**  radius  (R  =  Rn)« 

The  N  and  M  loops  in  the  computer  program  are  correspondingly 
named,  and  carry  out  the  determination  of  these  incremental  probabil¬ 
ities,  from  the  basic  relationship 

AP(N,  M)  »  P^  PdET  ^SL  ^L  ^POL 
where: 

P^  s  probability  of  downed  object  landing  in  the  area. 

Pqet  *  probability  of  detection  by  search  vehicle,  if  in  area 

searched. 

Pg^  ■  probability  of  safe  landing  of  downed  object. 

P|_^  >  probability  of  survival  (living  until  retrieved)  if  landing 

was  safe. 
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^POL  *  political  probability  of  retrieval  being  permitted. 
Further  let: 


t 


Tscr,  S 
■^SCR,  R 
"^END,  S 
"^END,  R 


=  time  after  impact,  running  variable. 

-  scramble  time  for  search  vehicle. 

=  scramble  time  for  retrieval  vehicle. 

=  endurance  time  of  search  vehicle. 

=  endurance  time  of  retrieval  vehicle. 

=  time  of  finding  (time  search  vehicle  will  reach  the  area). 
=  total  time  to  retrieval  for  that  area. 

=  speed  of  search  vehicle. 

=  speed  of  retrieval  vehicle. 

=  width  searched  on  either  side  of  vehicle's  path,  i.  e. ,  the 
search  half-width. 


^Ll  ”  viability  time  for  downed  airman  (or  object)  up  to  time 

of  finding. 

'  viability  time  from  finding  until  retrieval. 

In  this  model,  we  consider  that  FdET  ^  direct  function  only  of  the 
search  base  from  which  the  search  vehicle  leaves,  i.  e. ,  we  specify  the 
parameters  of  a  vehicle  through  its  base,  of  which  we  consider  it  to  form 
a  part.  The  same  is  true  of  W. 


The  probabilities  PsL  ^POL  direct  functions  of  the  territory 
where  the  object  is  downed,  i.  e. ,  of  the  global  sector  involved. 

Py^  is  a  direct  function  of  the  statistical  location  of  the  object  with 
respect  to  the  MPI,  and  also  indirectly  of  the  parameter  W  (search  half¬ 
width)  which  is  a  function  of  the  search  base  (since  this  includes  the 
search  vehicle,  as  noted  above). 


The  remaining  probability,  P^^,  is  a  direct  function  of  the  parameters 
Tj^j  and  ^1^2  ^^ich  are  characteristics  of  the  global  sector.  P]_^  is  also  a 
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function  of  the  times  of  finding  and  retrieval  since,  assuming  as  conven* 
ient  and  plausible  that  the  probability  of  remaining  alive  (the  viability  of 
the  object)  falls  off  exponentially  with  time,  we  have 

-Tp/T^i  -(Tf-Tt)/Tl2 

Pl  ■  «  X  e 

The  calculation  of  Tp  and  Tj  is  an  essential  problem  and  will  re¬ 
quire  all  the  remaining  parameters,  as  well  as  several  of  those  also 
entering  into  the  determination  of  other  probabilities. 

Note  that  Pqpj  has  been  specified  as  a  function  of  search  base 
(vehicle)  alone,  although  it  depends,  as  does  W,  on  the  characteristics 
of  the  downed  object.  To  attempt  to  include  explicit  dependence  on  the 
object  in  the  model  appears  unnecessarily  complicated  for  the  present 
purpose. 

It  should  also  be  noted  that  the  possible  delays  in  starting  search  or 
retrieval  have  been  considered  as  the  respective  scramble  times,  and 
that  the  endurance  times  of  the  vehicles  have  been  specified  directly  in¬ 
stead  of  as  ranges. 

Since  the  search  and  recovery  model  is,  by  definition,  global  in 
scope,  the  calculation  of  positions,  distances,  angles,  etc, ,  must  be  by 
a  method  suitable  for  such  a  problem.  Therefore,  we  consider  the  earth 
to  be  spherical  and  specify  all  positions  by  latitude  and  longitude  and  all 
distances  in  nautical  miles.  All  paths  of  vehicles  are  taken  as  along 
great  circles,  and  the  methods  of  spherical  trigonometry  were  originally 
employed,  A  nautical  mile  is  taken  to  be  the  length  of  one  minute  of  arc 
along  the  equator. 


DISTANCE  FORMULA 

For  the  distance  (Dist)  between  two  points  we  use,  in  most  cases, 
the  cosine  law  for  sides  from  spherical  trigonometry: 

.1 

Dist  ■  c  «  cos  (cos  a  cos  b  4  sin  a  sin  b  sin  C), 

We  always  take  this  distance  as  positive,  and  between  0  and  ir  radians 
inclusive,  so  that  no  ambiguity  results.  We  use  this  formula,  for  ex¬ 
ample,  in  choosing  bases,  since  any  of  the  criteria  call  for  comparing 

*lt  aotumlly  pnvdd  store  d0air«bl9  in  some  situations  to  miploy  the  methoda  ot  vector  algebra  to 
amid  dUtia^tiea  with  speoiaf  cases. 
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either  the  distance  from  the  MPI  to  each  base  or  other  quantities  de¬ 
pendent  on  this  distance. 

In  locating  the  point  F  (the  position  of  the  retrieval  vehicle  at  time 
Tp),  a  solution  by  spherical  trigonometry  is  possible,  but  a  solution 
using  vector  algebra  is  more  useful  since  it  avoids  numerous  ambiguities 
of  the  spherical  trigonometry  solution. 

SOLUTION  BY  SPHERICAL  TRIGONOMETRY.  Referring  to  Fig.  1, 
we  solve  as  follows: 


X-MPI 


F  -  Position  of  Rotriovol  Vohido  at  Tf 
B  -  Rotriovol  Bose 


P  -  I'lorth  Polo 


Fig.  1:  Spherical  Triangles,  Showing  Notations. 
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First,  in  large  -triangle  BXP,  we  have  given: 

a^  =  colatitude  (colat)  of  B,  i*  e. ,  of  retrieval  base 
cj  =  distance  from  retrieval  base  to  MPI 
b  =  colatitude  of  X,  i.  e« ,  of  MPI 
^  =  longitude  (long)  of  B  -  longitude  of  X 

and  we  want 

=  /-BXP. 

We  first  put 

®1  "  2  (aj+b+Cj). 

Then  from  the  tan  of  half-angle  law, 

^A  a  2  tan 

Second,  in  smaller  triangle  FXP,  we  have  given: 

^A 

'  ®1  ■  (Jf‘'^SCR,  r} 

and  we  want 

^€2  *  A  long  of  F 
a2  =  colat  of  F. 

Again  defining 

S2  =  i  (a2+b+C2) 
we  have 

^€2*2  tan"^ 

with  sign  of  /-C^. 


/8in(8i-b)  sin(8i-c71 
Ij  sinCsj-a)  sin  8^ 
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Then 

colatitude  of  F  >  a2 

longitude  of  F  s  (longitude  of  X4  Mod  2it. 

This  equation  holds,  however,  only  when  none  of  the  terms  in  the 
denominator  vanish.  These  cases  require  special  treatment.  There  are 
so  many  possibilities  that  the  program  becomes  inelegant  -  it  is  too 
complicated  for  proper  checkout  of  each  case.  Clearly  this  method  in¬ 
troduces  artificial  singularities  which  need  not  occur. 

SOLUTION  BY  VECTOR  ALGEBRA.  In  place  of  the  previous  nota¬ 
tion,  we  shall  now  define  the  unit  vectors  B,  f",  whose  origins  are 
at  the  center  of  the  earth  and  whose  termini  define  the  points  X,  B,  F, 
and  N  (North  Pole)  on  the  unit  sphere  (see  Fig.  2).  We  now  know  that 
vector  f'  lies  in  the  plane  of  ^  and  S,  thus  in  general,  it  is  true  that 


F  =  o  X  B  +  Y  (X  X  B). 

But,  since  in  this  problem,  X,  B,  and  F  are  coplanar, 

F  =  a  X  B 

U  ,  there  is  no  term  in  (X  k  B)* 

Note  that  we  shall  denote  the  "cross"  or  vector  product  as  (X  x  B)  and 
the  "dot"  or  scalar  product  by  (X*B). 

We  also  know 

F-B  =  cos  [vRTgj] 

F*  X  =  cos  [Dist  (B  to  X)  -  cos  (Vj^T^^)] 

B  'X  s  cos  ^Dist  (B  to  X)j. 

Then,  on  forming  these  products  and  solving  for  the  scalar  coefficients 
a  and/j,  we  find 

^  -F»X-  fB»jb(F.  B) 

D 

>  F»B  .  (F«X)(B«X) 
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where 

D  *  1  -  (B*X)2. 

Then,  in  the  general  case 

colat  F  =  cos”^  ja  cos  (colat  X)  +  cos  (colat  B)J 

and 

long  F  =  tan“^  [a  sin  (colat  X)  sin  (long  X) 

-t-yS  sin  (colat  B)  cos  (long  B)]. 


H  =  North  Polo 


Fig.  2:  Unit  Vectors  Used  in  Vector  Algebra  Solution, 
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The  only  situations  in  which  this  general  formula  cannot  be  used, 
are  those  where  D  =  0,  i,  e, ,  where  5*5^  =  ±1,  This  can  only  occur  for 
X  3  and  X  =  -B,  i,  e. ,  when  the  points  X  and  B  coincide,  and  when 
they  are  at  opposite  diameters  of  the  earth,  respectively.  The  first 
case  is  trivial;  the  second  allows  f'  to  be  on  any  great  circle  through  X 
and  B,  or  any  plane  through  X  and  S  as  we  express  it  in  this  solution. 

F or  our  problem,  we  may  pick  a  plane  arbitrarily  and  do  so  by  imposing 
some  further  condition.  Therefore,  in  this  case,  we  assume  that  the 
path  passes  over  the  North  or  South  Pole.  This  then  also  determines  a 
rather  special  situation,  and  the  appropriate  formulas  are  easily  written. 


SEARCH  CALCULATIONS 

The  basic  search  pattern  which  was  originally  contemplated  was  a 
spiral  of  Archimedes,  with  turns  spaced  by  2W,  i.  e. , 


where  R  =  radius  from  MFI,  6  s  angle  from  north  clockwise,  and  W 

3  sweep  half»width. 


Fig.  3:  Search  Pattern  Used  in  Model. 
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In  actually  constructing  a  model  to  be  mechanized,  it  was  found 
more  convenient  to  assume  a  search  by  concentric  circles  (see  Fig.  3), 
since  the  area  searched  in  a  given  time  is  essentially  the  same  and  the 
computations  are  simpler. 


We  assume,  therefore,  that  the  search  vehicle  proceeds  from  its 
base  directly  toward  the  MFI  of  the  downed  object,  but  upon  approach- 
ing  to  within  a  distance  W,  turns  and  makes  a  circle  about  the  MPI  at 
a  radius  of  W,  then  goes  out  radially  to  the  next  radius  which  is  2W 
greater,  makes  that  circle,  etc.  On  each  circle,  we  choose  eight 
points  at  angles  of  7r/4  to  each  other,  each  serving  as  a-representative 
sample  of  the  global  characteristics  in  its  general  area.  Therefore, 
we  find  its  coordinates  and  then  extract  the  parameters  pertinent  to 
our  problem  to  calculate  the  incremental  probability  for  that  area. 


COORDINATES  OF  CENTERS  (A)  OF  AREAS  SEARCHED.  To  cal¬ 
culate  the  coordinates  of  the  centers  (A)  of  the  areas  searched,  we  find 
points  about  the  MFI,  on  successive  radii  Rj^,  at  angles  0]^  from  north. 


Thus 

Rn  *  1  4  2W 


where 

Rj  =  W  or  Rq  =  -W 


and 


®N  =  ®N- 1  - T 
where 


01  =  -  Jy  or  00  =^. 

This  value  of  0^  was  chosen  so  that  R^^  will  never  pass  through  the 
pole,  unless  the  MPI  is  at  a  pole,  which  situation  we  treat  separately. 
In  general  (see  Fig.  4),  if  given 

0,  R,  and  a. 
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we  find 

c  =  cos  "  ^  ^cos  a  cos  r  +  sin  a  sin  r  cos  vj  , 

Then  defining 

s  *  i  (a+r+c) 
we  have 

A  L  =  2  tan  “  *  /sinfs-al  8in(s.-cJ 
V  sin  (s-r)  sin  s 

with  the  sign  of  (ir-6). 

The  longitude  A  =  longitude  X  +  AL  colatitude  A  =  c. 

The  only  special  cases  are  for  X  (the  object)  at  the  pole,  and  these 
are  treated  appropriately. 

The  area  probability  is  defined  as  the  probability  of  the  downed  ob¬ 
ject  lying  within  one  of  the  eight  sectors  of  the  ring  bounded  for  R  =  Rn 
by  the  concentric  circles  of  radii  Rj^+W  and  R^^-W, 

We  find  the  probability  of  the  object  lying  in  this  ring  as  the  proba¬ 
bility  of  lying  within  the  outer  disk  of  radius  (Rj^+W)  minus  the  probabil¬ 
ity  of  lying  within  the  inner  disk  of  radius  (Rj^-W),  Thus  defining  the 
outer  disk  probability  (Prob)  as 

2 


i/Rn+wV 


Disk  Probj 
and  the  inner  disk  probability  as 


Inner  Disk  Prob 
we  have 


Ring  Probj,j  =  Disk  Prob^  -  Inner  Disk  Prob^,. 


N 


N* 


But,  by  our  iteration  law 

Rj^  -  W  =  Rj^.  j  +  2  W  -  W  *  Rn.  1  +  W 
whence 


Inner  Disk  Probfj  *  Disk  Prob 


N-1 
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M  =  i'lorth  Pole 


Fig,  4:  Spherical  Triangle  Used  for  Finding  Point  A. 


and  thus 

Ring  Probj^  =  Disk  Probjq  -  Disk  Prob  j 
and  finally 

»  „  ^  Ring  Frobi^  _  Disk  Probvr  -  Disk  ProbM  i 

Area  Probj^  =  - - **  -  - 

We  also  note  that,  to  start  the  process  of  iteration,  we  have 
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Disk  ProbQ  =  0 


also 


<r  =  0.8493  x  CEP 


and  the  value  of  the  CEP  is  given  in  the  object  table. 

TIME  OF  FINDING.  The  time  of  arrival  of  the  search  vehicle  at  a 
given  area,  or  the  time  of  finding,  is  taken  as  a  constant  for  all  areas 
on  a  given  ring,  specifically  as  the  time  halfway  around,  i.  e. ,  as 


T  =  T 
F,  N  ^  F,  N-  1 


^  ^(^N  ^N- 1^ 


2W 


and,  for  the  first  ring 

Dist  (base  to  MPI)  -  W+irW 
^F,  1  ■  “^SCR,  S  V 

s 


Thus,  by  subtraction 


Dist  +  (it-3)W 
'^F,  0  =  "^SCR,  S  y 

s 

which  is  used  for  programming  convenience  as  to  starting  statement. 

ENDURANCE  TIME.  Before  these  area  probabilities  are  calculated, 
a  test  is  made  to  determine  whether  the  endurance  time  is  sufficient  for 
traversing  this  circle  of  radius  Rj^  by  testing  if 

Dist  +  ir  R  -  R 

"^F  ■  ^SCR,  S  ^  ^  ^END,  S  * 

8 

That  is,  can  the  circle  be  completed  and  can  return  to  base  be  made 
within  the  search  vehicle  endurance  time  ? 

TOTAL  AREA  SEARCHED.  The  total  area  searched  is  computed 
from  the  formula 


Rkt  +  W 
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where 

Ag  =  area  searched 
Rj^  =  radius  of  circle 
Rg  =  radius  of  the  earth. 


SECTION  4  -  DESCRIPTION  OF  COMPUTER  PROGRAM 


The  nature  of  the  model  which  is  mechanized  has  been  described 
earlier  in  this  report  and  the  formulas  used  have  been  given  in  detail  in 
Section  3  -  Mathematics  of  The  Model.  In  the  interest  of  clarity,  the 
flow  charts  in  the  following  description  deviate  from  the  actual  details  of 
the  program  in  a  few  places,  not  however,  to  an  extent  that  should  cause 
any  confusion  about  the  flow  of  the  program. 

The  APGC  modification  of  the  "FIB"  monitor  was  employed.  This 
permits  the  intermingling  of  main  programs  and  subprograms,  which 
may  individually  be  written  in  either  FORTRAN  or  FAP,  This  flexibility 
of  language  proved  most  advantageous  in  constructing  the  program.  The 
main  program,  including  tests  for  recognizing  types  of  input  cards,  was 
written  in  FORTRAN,  as  was  the  whole  of  the  computation  block.  How¬ 
ever,  several  subroutines  written  in  FAP  were  employed  where  packing 
and  unpacking  or  other  manipulation  of  bits  within  a  word  was  required. 


SIZE  AND  SPEED  OF  COMPUTER  PROGRAM 

The  program  proper,  including  necessary  library  subordinates,  occu¬ 
pies  roughly  8000  addresses  in  core,  and  the  tables  occupy  roughly 
another  8000,  so  that  in  all  about  16,  000  of  the  32,  768  (32  K)  addresses 
in  the  APGC  IBM  7090  installation  are  utilized  in  running  this  program. 

Of  the  addresses  taken  up  by  tables,  approximately  6000  are  used  for  the 
one  table,  GLOBE,  which  is  thus  a  major  item. 

It  takes  roughly  two  minutes  to  compile  the  program  from  the  sym¬ 
bolic  decks.  The  time  required  to  actually  compute  a  factorial  experi¬ 
ment  will  naturally  vary  considerably,  but  in  preliminary  testing,  32 
case  experiments  which  appeared  typical  were  run  in  about  1  1/2  minutes. 
Printing  was  done  off-line  using  an  IBM  1401, 
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Flow  Chart  for  Overall  Frogram, 
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DUPLICATIOJSr  OF  CERTAIN  TABLES 

In  order  to  carry  out  the  variation  of  parameters  called  for  by  a 
factorial  experiment  it  is  necessary  to  vary  the  factors  independently. 

In  fact,  it  is  necessary  to  restore  each  table  to  its  original  or  standard 
values  as  set  up  before  applying  each  modification,  (Certain  "perma¬ 
nent"  tables  are  not  subject  to  such  variations,)  It  was  thus  necessary  to 
duplicate  those  tables  which  can  be  factorially  modified;  one  set  (the 
upper  tables)  contained  the  SETUP  or  standard  values  for  resetting  pur¬ 
poses,  the  other  set  (the  lower  tables)  contained  the  current  values  to  be 
used  with  those  in  the  permanent  tables  in  calculating  a  given  case. 

In  actual  program,  all  tables  except  the  permanent  tables  (GLOBE, 
VIAB,  and  PROB)  are  duplicated.  Thus,  in  addition  to  the  tables  SBASE, 
RBASE,  OBJ,  SDOCT,  and  RDOCT,  which  will  now  be  referred  to  as 
lower  tables,  there  is  also  another  set  of  tables,  of  identical  format, 
which  are  known  as  upper  tables.  Within  the  program,  these  upper 
tables  are  called  USBASE,  URBASE,  etc, ,  though  not  so  prefixed  by  U 
in  input  card  formats  for  either  SETUP  or  FACTOR  cards.  See  Appen¬ 
dix  II  for  a  discussion  of  input  card  formats. 


MAIN  PROGRAM 

The  heart  of  the  program  is  block  1000,  the  computation  block. 

Every  other  part  deals  with  manipulations  of  tables  and  lists,  setting  up 
various  combinations  of  conditions  under  which  the  computation  block  will 
be  executed.  Fig,  5  is  the  flow  chart  of  the  overall  program;  Fig.  6  is  a 
more  detailed  chart  of  the  computation  block.  Except  where  noted,  the 
program  was  written  in  FORTRAN.  The  numbers  of  the  following  para¬ 
graphs  correspond  to  the  appropriate  block  in  Fig.  5,  and  the  function  of 
each  block  is  described. 

100,  Certain  often-used  constants  are  given  names  and  page  head¬ 
ings,  and  numbering  blocks  are  initialized. 

105.  An  input  card  is  read.  It  is  later  interpreted  in  block  109.  If 
no  card  is  present  when  control  reaches  block  105,  control  returns  to 
the  monitor  and  the  machine  run  is  over, 

108.  Each  input  card  is  immediately  printed  to  provide  a  more  com¬ 
plete  reconstruction  of  the  computation  for  the  output. 

120,  125,  109,  1  10,  130,  115.  A  seri'»s  of  tests  is  made  to  determine 
which  type  card  was  just  read  in,  and  control  is  sent  to  the  proper  block, 
i.  e.  ,  to: 
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300  -  if  a  LOAD  card  is  being  read. 

400  -  if  an  UNLOAD  card  is  being  read. 

600  -  if  a  GOGO  card  is  being  read. 

500  -  if  a  SETUP  card  is  being  read. 

200  -  if  a  FACTOR  card  is  being  read. 

800  -  if  a  DUMP  card  is  being  read. 

135  -  if  card  format  is  not  one  of  the  above,  i.  e. ,  incorrect. 


200.  A  test  is  made,  by  a  FAP  subroutine  (SUBl)  to  determine  what 
sort  of  FACTOR  card  is  present,  and: 


a.  If  word  2  (cols  7.12  incl)  consists  of  RESET,  control  passes 

to  250. 


b«  If  word  2  has  allowable  factors  and  levels  (shown  in  cols  8 
and  10  through  12  respectively),  control  passes  to  205. 

c.  If  word  2  does  not  meet  either  criteria,  the  card  is  ignored, 
and  control  passes  to  105  and  another  card  is  read. 

205.  The  information  is  stored  in  the  FACTOR  list  in  the  position 
designated  by  the  indices  L  and  N,  M. 

250.  The  FACTOR  list  is  reset.  A  NULL  factor  is  stored  as  level 
No.  1  of  each  factor  and  DONE  is  inserted  for  all  other  levels  including 
the  1 1^^  levels,  which  can  contain  no  other  data. 

300.  The  tables  are  read  into  memory  from  magnetic  tape.  The 
tables  must  have  been  generated  and  stored  on  previous  runs  of  this  pro¬ 
gram  by  the  use  of  SETUP  cards  and  an  UNLOAD  card.  (See  blocks  400 
and  500. ) 

400.  The  tables  are  written  on  magnetic  tape  for  the  purpose  of 
being  used  on  subsequent  runs  of  this  program,  as  described  under  block 
300  above. 

500.  The  SETUP  cards  are  used  to  insert  values  into  the  various 
upper  and  permanent  tables.  A  FAP  subroutine  (MOD)*  is  used  to 


*MOD  Thia  is  a  mibroutioe  writtm  in  FAP  wkidi  im  vofioyad  to  "modUy**  my  table  in  memory. 
The  ^‘Modification"  includea  original  netting  up  of  the  data  in  theae  tablea.  It  ia  called  up  by 
either  block  500  or  block  700. 
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decode  the  necessary  information  and  to  decide  which  table  and  which 
portions  of  it  are  to  be  set  up. 

600.  This  block  and  block  605  are  interconnected  in  the  program. 
Together  they  pick  which  of  the  factorial  cases  is  to  be  run  next.  At 
block  600,  the  process  is  started  by  choosing  the  first  level  for  each  of 
the  five  factors,  and  transferring  control  to  block  700.  After  the  first 
case  has  been  computed  by  block  1000,  control  returns  to  605. 

605.  We  now  wish  to  choose  the  next  case.  We  do  this  by  stepping 
up  factor  5,  and  if  the  next  level  is  not  a  DONE,  transferring  control  to 
block  700,  If  this  higher  level  is  a  DONE,  however,  this  factor  is  re¬ 
stored  to  its  first  level,  and  the  next  lower  factor  is  stepped  up  and 
again  tested.  This  is  continued  until  either  some  factor  is  stepped  into 
a  level  not  containing  a  DONE,  with  resulting  transfer  of  control  to  700, 
or  it  is  determined  that  the  entire  set  of  factorial  combinations  has  been 
run  through,  whereupon  control  passes  to  105. 

700.  This  block  carries  out  the  final  readying  for  actual  computa¬ 
tion  by  ordering  the  modifications  in  the  lower  tables  which  correspond 
to  the  information  in  the  factor  and  level  chosen  by  block  605.  Actually 
this  modification  is  carried  out  by  subroutine  MOD. 

Thus  actually  for  each  case,  block  700  first  restores  all  lower 
tables,  then  makes  the  modifications  required  by  block  605  (or  600), 
operating  through  MOD,  as  noted. 

When  called  by  block  500,  MOD  is  used  to  make  the  modifications 
in  permanent  tables  or  in  upper  tables  as  called  for  by  SETUP  cards. 

When  called  by  block  700,  MOD  is  used  to  make  modifications 
to  lower  tables  as  directed  by  that  block,  i.  e. ,  it  resets  all  lower  tables 
to  the  values  in  their  corresponding  upper  tables,  then  modifies  them  as 
required  by  the  particular  factor  table  combination  chosen  in  block  605 
(or  600), 


COMPUTATION  BLOCK 

The  detailed  flow  chart  of  the  computation  block  is  given  in  Fig,  6, 
This  block  is  written  entirely  in  FORTRAN,  Control  reaches  this  block 
from  700  after  all  specified  factorial  modifications  have  been  made.  The 
computer  block  consists  of  three  portions.  The  first  deals  with  certain 
preliminary  calculations  including  the  selection  of  the  search  vehicle  and 
of  the  retrieval  vehicle.  The  second  deals  with  the  generation  of  the 
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representative  points,  and  the  calculation  of  the  quantities  associated  with 
them.  The  third  computes  some  additional  values  and  prints  the  one  line 
of  output  generated  by  the  whole  process. 

The  second  portion  of  the  computer  block  consists  of  two  nested 
loops.  The  outer  (N)  loop  is  associated  with  the  radii  of  the  series  of 
concentric  circles  about  the  MPI,  The  inner  (M)  loop  is  associated  with 
the  angular  location  of  the  series  of  representative  points  lying  on  the 
circumferences  of  these  circles.  The  M  loop  is  executed  exactly  eight 
times  (eight  equally  spaced  points  are  generated  on  each  circle),  and 
the  N  loop  can  be  executed  a  maximum  of  200  times,  but  is  almost  always 
left  earlier  due  to  some  condition  of  time  or  area  covered, 

1000.  Certain  accumulated  quantities  are  set  to  zero,  and  infor¬ 
mation  about  the  probable  location  of  the  object  and  the  method  of  search 
and  recovery  to  be  used  is  extracted  from  tables  OBJ,  SDOCT,  and 
RDOCT,  respectively.  This  consists  of  the  colatitude  and  longitude  of 
the  MPI  and  two  numbers  denoting  the  criteria  for  choosing  the  search 
base  and  retrieval  base. 

1002.  The  best  search  base  is  chosen  from  the  table  SBASE  accord¬ 
ing  to  the  search  doctrine  DRITS.  DRITS  determines  whether  the  chosen 
base  is  to  be  best  in  the  terms  of;  (1)  closest  to  the  MPI,  (2)  can  reach 
the  MPI  soonest,  (3)  having  reached  the  MPI  can  search  longest  or,  (4) 
can  search  the  greatest  area. 

1011.  It  is  determined  here  whether  the  best  search  vehicle  has  the 
endurance  to  reach  the  MPI  and  have  at  least  some  time  left  for  search¬ 
ing  before  it  is  necessary  to  return.  If  not,  the  calculation  is  terminated. 

1015.  It  is  determined  whether  retrieval  is  to  be  accomplished  by 
the  search  vehicle  itself  or  by  a  separate  retrieval  vehicle.  This  is 
under  the  control  of  KRITR  (retrieval  criterion)  obtained  from  RDOCT 
in  block  1000. 

10  16.  A  retrieval  base  is  chosen. 

1024.  It  is  determined  whether  to  terminate  the  calculation  due  to 
the  lack  of  an  adequate  retrieval  vehicle  (one  that  can  make  it  to  the 
MPI  and  back  to  its  base).  Information  is  determined  as  to  the  locations 
of  the  retrieval  base,  its  distance  from  the  MPI,  and  the  speed,  endur¬ 
ance  time,  and  scramble  time  of  the  retrieval  vehicle. 

1030.  Here  the  initialization  for  the  N  (ring)  loop  is  performed, 

Rp^,  the  radius  of  the  circle,  is  set  so  that  the  first  circle  will  have 
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radius  W,  the  sweep  width  of  the  search  vehicle,  and  Tp  is  set  so  that 
its  value  for  the  first  circle  will  be  the  average  time  of  finding  the  object 
if  it  is  in  the  first  ring. 

1031.  Here  R  is  increased  by  2W  so  that  successive  circles  will  be 
spaced  in  such  a  way  that  the  search  vehicle  can  cover  all  the  area  be¬ 
tween  them.  Tp  is  increased  by  the  appropriate  amount  so  as  to  repre¬ 
sent  the  average  time  it  would  take  to  find  the  object  if  it  is  in  the  current 
ring. 


1032.  It  is  determined  whether  sufficient  time  remains  for  the  search 
vehicle  to  complete  the  search  of  the  current  ring  and  return  to  its  base 
within  its  endurance  time.  If  not,  search  is  terminated. 

1034.  The  probability  that  the  object  is  in  one  of  the  eight  sectors 
of  this  ring  is  calculated.  This  probability,  or  AREAP,  will  be  used 
as  a  weight  in  the  weighted  average  of  individual  probabilities,  and  is 
calculated  at  this  point  of  the  program  since  it  is  the  same  for  each  sec¬ 
tor  of  this  ring. 

1036.  If  a  separate  vehicle  is  used,  we  next  calculate  its  expected 
position  at  the  time  of  finding  the  object  if  the  object  is  found  in  this  ring. 
This  point  (F)  is  determined  from  Tp,  the  position  of  the  MPI,  and  infor¬ 
mation  about  the  location  of  the  retrieval  base  and  the  capabilities  of  the 
retrieval  vehicle  obtained  in  block  1016.  It  is  assumed  that  when  the 
search  vehicle  sights  the  object,  it  will  inform  the  retrieval  vehicle  of 
the  object's  location  so  that  the  retrieval  vehicle's  course  can  be  altered 
to  proceed  directly  to  the  location.  This  point  of  turning  is  needed  later 
to  calculate  the  time  the  retrieval  vehicle  reaches  the  object.  This 
point  is  assumed  to  be  the  same  for  each  sector  of.  this  ring.  If  no  sepa- 
i  -vte  retrieval  vehicle  is  used,  this  portion  of  the  program  is  bypassed. 

1065.  Here,  initialization  takes  place  for  the  M  (or  angle)  loop. 

9  starts  from  0,  2  radian  rather  than  from  0  so  that  A  is  never  due  north 
of  the  MPI.  This  has  the  effect  of  eliminating  certain  exceptional  condi¬ 
tions  from  the  subsequent  calculation  of  the  position  of  A. 

1066.  9  is  stepped  by  ir/4,  and  the  coordinates  of  A  are  calculated. 
Thus,  the  eight  points  will  be  equally  spaced  about  the  circle. 

1075.  In  this  block,  the  position  of  A  is  used  to  determine  which 
cell  in  the  GLOBE  table  contains  the  description  of  the  area  containing  A, 
and  this  information  is  unpacked. 

1080.  At  this  point  all  the  necessary  information  has  been  obtained 
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to  calculate  the  probability  of  successfully  retrieving  the  object  if  it  is 
at  the  point  A.  This  is  done,  after  which  it  is  multiplied  by  the  weight¬ 
ing  factor  to  find  the  incremental  probability  AF(N,  M),  and  accu¬ 
mulated  to  give  the  cumulative  probability  F  for  the  points  so  far  tested, 

1 130.  If  all  eight  points  A  have  been  generated,  the  program  pro¬ 
ceeds.  Otherwise  control  is  returned  to  the  top  of  the  M  loop  to  generate 
another  angle. 

1 132.  At  this  point,  it  is  determined  whether  the  probability  that 
the  object  is  inside  the  area  searched  is  greater  than  0.  999  on  the  basis 
of  the  CEP  and  the  current  value  of  R|^.  If  so,  the  generation  of  rings 
is  terminated. 

1 134.  If  control  reaches  this  point  and  less  than  200  rings  have  been 
generated,  control  is  returned  to  the  top  of  the  N  loop  to  generate  another 
radius.  Otherwise,  ring  generation  is  terminated. 

1 1 38.  KTERM  is  a  flag,  denoting  the  reason  the  generation  of 
points  A  was  terminated.  This  information  is  printed  as  output  for  the 
benefit  of  the  user.  For  this  program,  it  is  set  to  2  to  denote  that  the 
maximum  allowable  number  of  areas  (8  x  200)  were  covered. 

1 140.  KTERM  is  set  to  1  to  denote  that  calculation  was  terminated 
due  to  reaching  the  end  of  the  search  vehicle's  endurance  time. 

1  142,  KTERM  is  set  to  3  to  denote  the  lack  of  a  suitable  search  ve¬ 
hicle,  i.  e. ,  one  which  could  reach  the  MPI  and  return  within  its  endur¬ 
ance  time. 

1 144.  KTERM  is  set  to  4  to  denote  the  lack  of  a  suitable  retrieval 
vehicle,  again  one  with  sufficient  endurance  to  reach  the  MPI  and  return 
to  base. 


1 150.  The  total  area  in  square  miles  is  calculated  from  the  final 
value  of  Rj^. 

1200.  At  this  point,  all  output  quantities  have  been  calculated  and 
are  printed.  A  case  has  now  been  completed,  and  control  passes  to 
block  605. 
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SECTION  5  -  DISCUSSION  OF  RESULTS 


As  has  been  described,  a  model  simulating  a  specific  type  of  Global 
Search  and  Recovery  problem  with  variable  parameters  was  actually  con¬ 
structed  and  tested.  Considerable  insight  was  gained  into  the  nature  and 
magnitude  of  the  problems  which  would  be  encountered  in  any  such  simu¬ 
lation  program.  The  global  requirement  was  the  determining  factor  in 
the  problem. 

The  concept  of  characterizing  the  surface  of  the  globe  by  dividing  it 
into  regions  within  which  the  several  pertinent  parameters  could  be  con¬ 
sidered  constant  still  seems  a  basically  valid  one.  However,  the  prac¬ 
tical  difficulties  of  choosing  a  satisfactory  size  for  a  region  (called  sec¬ 
tor  in  this  report)  and  of  actually  coding  such  a  global  table  for  a 
realistic  problem  were  seen  to  be  enormous.  Even  the  coding  of  the 
relatively  small  number  of  sectors  chosen  for  this  investigation  was 
soon  found  to  be  far  from  a  simple  chore.  Such  difficulties  appear  in¬ 
herent  to  the  problem,  and  no  real  solution  is  evident.  Probably  spe¬ 
cific  choices  of  area  size  will  be  required  depending  on  the  nature  of 
the  simulation  to  be  perfornied  and  on  the  available  resources  for  ascer¬ 
taining  the  pertinent  parameters  and  for  actually  coding  the  global  sector 
table  to  be  employed.  For  limited  areas,  the  problem  can  of  course  be 
greatly  simplified,  and  this  may  be  required  for  nonglobal  problems. 

The  concept  of  assigning  the  remaining  parameters  to  search  or 
retrieval  bases  or  to  the  searched-for  object  itself  was  also  a  simpli¬ 
fication  of  actual  conditions.  However,  it  seems  to  be  a  reasonable  con¬ 
cept,  at  least  for  preliminary  studies,  using  a  model  similar  to  that 
here  considered. 

The  model  itself,  with  its  single  type  of  search  and  two  fixed  types 
of  retrieval,  was  intended  primarily  as  an  example  of  the  sort  of  models 
which  might  be  employed.  It  is  too  inflexible  to  be  entirely  represent¬ 
ative  of  the  types  of  search  and  retrieval  problems  to  be  anticipated. 

In  future  work,  a  more  general  model,  utilizing  Monte-Carlo  techniques, 
appears  to  be  a  more  promising  approach.  Nevertheless,  the  model  used 
did  reveal  considerable  information  about  the  type  of  conceptual  and  even 
geometrical  problems  which  occur,  for  example,  the  difficulties  encoun¬ 
tered  with  ambiguities  and  special  cases  if  the  method  of  spherical  trig¬ 
onometry  are  used,  together  with  ways  of  avoiding  most  of  these  troubles. 
Also,  it  soon  became  clear  that  any  realistic  model  of  retrieval  using 
land  vehicles,  taking  account  of  constraints  such  as  road  systems  avail¬ 
able  and  topographical  barriers  to  nonroad  vehicle  travel,  could  easily 
become  extremely  complex.  For  limited  areas,  again  the  situation 
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might  be  quite  different.  Poss'bly  on  the  global  scale,  some  sort  of 
average  travelability  parameters,  analogous  to  those  for  probability  of 
safe  landing,  etc. ,  might  be  used,  again  in  broadbrush  fashion.  This 
remains  an  unexplored  problem. 

The  concepts  of  using  tables  for  parameters  and  of  the  factorial 
variability  of  these  parameters  appear  to  be  basically  sound  and  useful. 
In  view  of  the  extremely  large  number  of  global  sectors  for  which  pa¬ 
rameters  need  definition,  and  the  small  but  not  insignificant  number  of 
parameters  required  for  describing  other  conditions,  the  need  for  a 
means  for  introducing  such  parameters  in  a  convenient  straightforward 
manner  is  compelling.  It  must  be  done  so  as  to  minimize  the  effort  re¬ 
quired  and  hold  errors  to  a  minimum. 

As  examples  of  minor  shortcomings  of  the  present  program  which 
proved  annoying  in  this  respect,  the  following  is  noted: 

1.  The  width  in  longitude  of  the  global  sectors  was  chosen  as 
4  degrees,  to  give  what  seemed  to  be  reasonable  balance  between  size 
and  number  of  sectors.  This  choice  was  unfortunate  since  most  maps 
are  graduated  in  multiples  of  5  degrees. 

2.  The  formats  of  input  cards  used  to  introduce  parameters 
were  very  similar,  but  the  few  deviations  from  truly  identical  format 
which  were  made  caused  errors. 

3.  The  omission  of  a  printout  of  the  "factor  table"  proved  a 
handicap,  since  after  a  succession  of  modifications,  its  contents  were 
not  obvious. 


SECTION  6  -  CONCLUSIONS 


Simulation  of  Global  Search  and  Recovery  problems  can  never  be 
easy  and  is  possible  only  by  greatly  simplifying  such  factors  as  the  geo¬ 
graphy  of  the  earth  and  its  effects.  By  making  such  simplifying  assump¬ 
tions,  models  of  considerable  promise  can  be  constructed  which  are 
amenable  to  use  on  high  speed  computers.  The  degree  to  which  such 
models  will  be  realistic  remains  to  be  determined. 
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APPENDIX  I 

FACTORIAL  EXPERIMENTS 


DEFINITION  OF  FACTORIAL  EXPERIMENT 

The  term  "factorial  experiment"  comes  from  the  field  of  Statistical 
Design  and  Analysis  of  Experiments.  It  applies  to  an  experiment  in 
which  the  effects  of  a  number  of  different  factors  are  investigated  simul¬ 
taneously.  This  is  done  by  setting  up  the  complete  set  of  cases  (or 
"treatments"  in  the  statistician's  nomenclature)  comprising  all  combi¬ 
nations  of  all  "levels"  investigated  of  each  factor.  Thus  we  might  choose, 
for  example: 

FACTOR  I  =  Sweep  Width 

LEVEL  1  =  0.  5  NM 
LEVEL  2  *  10  NM 
LEVEL  3  =  100  NM 

FACTOR  II  =  CEP 


LEVEL  1  a  10  NM 
LEVEL  2  =  100  NM 

FACTOR  III  =  Retrieval  Doctrine 


LEVEL  1  =  Search  vehicle  also  retrieves 
LEVEL  2  =  Separate  vehicles  used  for  retrieval 


This  would  provide  a  possibility  for  12  (3  x  2  x  2)  cases,  as  follows: 


Case  Factor  I 
No.  at  Level 


Factor  II  Factor  III 
at  Level  at  Level 


1  1 

2  1 

3  1 

4  1 

5  2 

6  2 

7  2 

8  2 

9  3 

10  3 

11  3 

12  3 


1 

1 

2 

2 

1 

1 

2 

2 

1 

1 

2 

2 


1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 
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In  an  actual  experiment  it  is  necessary  to  repeat  (replicate)  each 
case  and  to  randomize  the  order  in  which  the  cases  are  run.  Neither 
of  these  requirements  holds,  naturally,  for  the  simulated  experiments 
described  in  this  report. 


POSSIBLE  METHODS  OF  VARYING  PARAMETERS 

ELEMENTARY  METHOD.  Factorial  experiments  could  be  run 
"bruteforce"  by  specLying  each  case  separately  using  SETUP  cards, 
followed  for  each  case  by  a  GOGO  card.  The  number  of  cards  required 
for  such  an  experiment  may  become  quite  considerable  for  sizeable 
experiments.  In  Table  I-l,  the  number  of  cards  so  required  is  com¬ 
pared  with  the  corresponding  number  required  using  the  Factor  Table 
method  described  below.  By  a  3  x  2  factorial  design,  for  example,  we 
mean  one  factor  at  each  of  3  levels  and  another  at  each  of  2.  This  is  a 
standard  notation. 

TABLE  I-l.  COMPARISON  OF  METHODS  OF  VARYING  PARAMETERS. 


Factorial 

Design 

Number  of  Modification  Cards  Required  | 

Number 
of  Cases 

Using  Additional 
SETUP  Cards* 

Using  FACTOR 
Cards 

III  ,111 

2x1 

2 

1 

1 

2x2 

4 

4 

2 

3x2 

6 

7 

3 

3x3 

9 

10 

4 

4x2 

8 

10 

6 

8 

11 

3 

27 

36 

6 

2x2x2x2 

16 

26 

4 

3x3x3x3 

81 

115 

8 

2x2x2x2x2 

32 

57 

5 

3x3x3x3x3 

243 

323 

10 

2x3x4x5 

120 

148 

10 

10x10x10 

1,000 

1,  107 

27 

10x10x10x10x10 

111,  105 

45 

*  In  addition  a  number  of  GOGO  cards  equal  to  the  number  of  cases 


would  be  required. 

*  In  addition  one  FACTOR  RESET  and  one  GOGO  card  are  required. 
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The  formula  for  the  number  of  cards  required  by  the  elementary 
method  is 


5 

Number  of  cards  =  M 

^1  ^ 


where  Mp  =  number  of  cards  to  modify  the  factor,  and  where 

Mp  =  U2*l2  *  *  *  ^  =  0  if  1  p  =  1,  and  Ip  =  number 

of  levels  in  P*^  factor. 

FACTOR  TABLE  METHOD,  In  this  method,  a  much  smaller  num¬ 
ber  of  cards  need  be  punched,  in  fact 

5 

number  of  cards  =  (Ip  -  1) 

P=1 


where  Ip  =  number  of  levels  of  the  P*^  factor  (only  one  card  for  each 
level,  except  the  first,  of  each  factor).  Here  we  have  assumed  through¬ 
out  for  comparison,  that  the  first  level  of  each  factor  will  be  the  un¬ 
changed  original  value,  which  is  the  normal  situation. 

It  should  be  noted  that  the  number  of  cases  which  will  be  computed, 
not  the  number  of  cards  required,  will  undoubtedly  be  the  limiting 
factor  in  running  factorial  experiments  by  this  method.  Even  for  a 
computer  such  as  the  IBM  7090,  an  experiment  with  100,  000  cases  will 
require  more  computing  time  than  is  likely  to  be  warranted. 

In  addition  to  the  tables  described  elsewhere,  the  program  also 
provides  a  region  in  memory,  FACT  (I,  J,  K),  which  is  essentially 
another  three-dimensional  table,  and  which  may  be  called  the  factor 
table.  Information  as  to  what  modifications  are  to  be  made,  factorially, 
in  the  main  tables  during  an  experiment  is  stored  in  this  region.  The 
indices  L  and  M  refer,  respectively,  to  the  Factor  and  to  the  level  of 
the  modification  for  which  data  is  stored.  Essentially  this  data  is  iden¬ 
tical  with  that  contained  in  a  SETUP  card,  except  that  it  cannot  refer  to 
GLOBE,  PROB,  or  VIAB.  It  is  used  preceeding  each  calculation  to 
modify  the  proper  lower  tables  after  they  have  been  restored  to  the 
standard  values  as  stored  in  the  upper  tables  by  the  SETUP  cards. 

The  third  index  (K)  is  required  only  because  all  the  data  which  can 
be  given  on  one  card  cannot  go  into  one  word  -  in  fact,  ten  words  are 
required.  For  simplicity,  this  index  may  be  neglected  and  the  factor 
table  as  being  a  two-dimensional  table  set  up  as  follows: 
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LEVEL 

FACTOR 

I 

II 

III 

IV 

V 

1 

FI,  1 

F2,  1 

F3,  1 

F4,  1 

F5,  1 

2 

FI,  2 

F2,  2 

F3,  2 

F4,  2 

F5,  2 

3 

4 

FI,  3 

F2,  3 

F3,  3 

F4,  3 

F5,  3 

• 

• 

10 

FI,  10 

F2,  10 

F3,  10 

F4,  10 

F5,  1( 

11 

DONE 

DONE 

DONE 

DONE 

DONE 

The  cue  to  the  program  that  all  specified  factorial  modifications 
have  been  made  for  a  given  factor  is  the  presence  of  the  word  DONE 
in  the  next  higher  level.  We  thus  provide  an  11*^  level  in  the  factor 
table  which  always  contain  DONE's  and  is  not  addressable  by  the  user 
of  the  program.  A  given  level  which  consists  of  no  change  is  indicated 
by  the  presence  of  the  word  NULL  in  that  level.  Each  level  of  each 
factor  (10  X  5  or  50)  can  be  set  by  a  separate  card. 

Factor  Reset  Card.  For  convenience,  we  have  provided  a  FACTOR 
RESET  card,  which  by  itself,  will  set  up  the  following  table: 


LEVEL 

FACTOR 

I 

II 

III 

IV 

V 

1 

NULL 

NULL 

NULL 

NULL 

NULL 

2 

DONE 

DONE 

DONE 

DONE 

DONE 

3 

• 

DONE 

DONE 

DONE 

DONE 

DONE 

• 

• 

10 

DONE 

DONE 

DONE 

DONE 

DONE 

11 

DONE 

DONE 

DONE 

DONE 

DONE 

This  is  the  desired  table  for  a  simple  single-case,  no- modification, 
factorial  experiment,  and  if  no  further  FACTOR  cards  precede  the  COCO 
card,  that  one  case  is  what  will  be  run  as  that  entire  experiment. 

General  Experiment.  The  FACTOR  RESET  card  is  also  useful  in 
general,  since  this  one  card  sets  in  the  required  DONE's  and  NULL's, 
and  only  those  modifications  desired  need  be  specified  by  other  factor 
cards. 


For  example,  suppose  only  the  following  factorial  variations 
are  desired  expressed  symbolically: 
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FACTOR  1, 

2 

Ml 

FACTOR  1, 

3 

Ml 

FACTOR  3. 

1 

M3 

2 

3 

1 


Then  the  table  configuration  resulting  from  use  of  FACTOR  RESET 


followed  by  the 

se  three 

cards  would  be, 

symbolically: 

LEVEL 

FACTOR 

I 

II 

III 

IV 

V 

1 

NULL 

NULL 

M3,  1 

NULL 

NULL 

2 

Ml,  2 

DONE 

DONE 

DONE 

DONE 

3 

M2,  3 

DONE 

DONE 

DONE 

DONE 

4 

• 

DONE 

DONE 

DONE 

DONE 

DONE 

• 

• 

10 

DONE 

DONE 

DONE 

DONE 

DONE 

11 

DONE 

DONE 

DONE 

DONE 

DONE 

Note  that  in  this  example  we  have  chosen  not  to  compute  a  case  v/ith 
Factor  III  at  its  setup  values  at  all.  This  can  be  done  although,  in 
general,  it  is  undesirable. 

The  11*^  level  DONE's  are  necessary  to  safeguard  the  program,  and 
higher  levels  have  no  meaning.  Any  attempt  to  store  factors  above  five, 
or  levels  above  10  will  be  ignored,  since  tests  are  provided  to  prevent 
such  assignments. 

What  is  chosen  as  a  given  factor  is  completely  flexible  except  for 
parameters  in  the  GLOBE  or  its  auxiliary  tables.  Different  levels  of  a 
factor  may  even  be  different  parameters.  However,  since  this  disturbs 
the  experiment  from  its  truly  factorial  character,  such  usage  is  not 
generally  recommended,  although  it  may  be  desirable  in  special  cases. 

In  actual  operation  when  a  GOGO  card  is  read  before  a  case  is  run, 
the  program  restores  all  lower  tables  to  the  upper  table  (standard)  values. 
It  then  makes  the  changes  called  for  by  that  case,  runs  it,  the  case  re¬ 
stores  the  valves,  then  modifies  lower  tables  again  until  all  cases  called 
for  have  been  run.  The  experiment  is  then  terminated,  and  the  program 
reads  another  card,  which  may  start  the  specification  of  the  parameters 
for  a  new  experiment  by  SETUP  and/or  FACTOR  cards  or,  if  no  further 
card  is  present,  terminate  the  machine  run. 
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APPENDIX  II 

INPUT  CARD  FORMATS  AND  AUXILIARY  TABLE  VALUES 


INPUT  CARD  FORMATS 

The  formats  of  the  various  input  cards  are  shown  in  Fig,  II- 1, 
Where  Ij  and  I^  and  J|  and  J2  are  shown  as  in  the  GLOBE  table,  all 
members  of  that  table  for  values  of  I  from  I]  to  I2  and  J|  to  inclusive 
are  to  be  affected.  For  SBASE  and  RBASE,  only  I^  and  I^  are  used  with 
a  similar  meaning.  Where  flags  are  shown,  the  inclusion  of  any  punch 
in  this  column  causes  entry  of  the  corresponding  identified  data  which 
otherwise  is  ignored.  For  FACTOR  cards,  F  and  LL  correspond  to 
factor  and  level  respectively. 


GLOBE  WORD  FORMAT 

As  a  matter  of  interest,  the  format  of  a  GLOBE  word  is  shown  in 
Fig,  II- 2, 

OUTPUT  FORMAT 

As  each  card  is  read  in,  it  is  printed  out  in  identical  format*  This 
will  appear  in  the  printed  output*  In  addition,  any  tables  for  which 
DUMP  cards  appear  will  be  printed  out  for  later  checking*  The  formats 
of  these  tables  will  not  be  described  here  but  can  be  found  in  the  original 
program* 

The  essential  output,  i*e*,  the  computed  results  from  each  experi¬ 
ment,  is  printed  out  after  the  printout  of  each  GOGO  card*  An  illus¬ 
trative  printout  for  a2x4x2x2xl,  i*e,,  32-case  experiment,  is 
shown  in  Fig*  II- 3*  This  shows  the  way  the  problem  was  set  up,  the 
,  factor  table,  and  the  computed  results  as  printed  out. 


SUGGESTED  VALUES  FOR  AUXILIARY  TABLES 

Suggested  values  for  the  PROB  auxiliary  tables  are  shown  in  Table 
II- 1,  and  the  values  for  VIAB  are  shown  in  Table  II-2*  The  subroutine 
FAPCON  is  required  as  it  is  used  by  MOD  in  decoding  the  values  in 
these  tables. 
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Fig,  II- 1:  Formats  of  Input  Cards 
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Fig,  II-2:  Globe  Word  Format, 


TABLE  II- 1,  PROBABILITY  (PR OB)  TABLE  (SUGGESTED  VALUES), 


K 

PROB 

xyz 

1 

,00 

000 

2 

.01 

3 

.02 

4 

.03 

5 

.04 

040 

6 

.05 

050 

7 

.07 

070 

8 

.08 

080 

9 

.  10 

100 

10 

.  12 

120 

11 

.  15 

150 

12 

.20 

200 

13 

.25 

250 

14 

.30 

300 

15 

.35 

350 

16 

.40 

400 

K 

PROB 

xyz 

17 

.45 

450 

18 

.50 

19 

.55 

20 

,60 

600 

21 

.65 

650 

22 

.70 

700 

23 

.75 

750 

24 

.80 

800 

25 

.85 

850 

26 

.90 

900 

27 

.95 

950 

28 

.96 

960 

29 

.97 

970 

30 

.98 

980 

31 

.99 

990 

32 

1.0 

101 
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TABLE  II-2,  VIABILITY  (VIAB)  TABLE  CODING  (SUGGESTED  VALUES), 


K 

xyz 

Time 

(Hr) 

K 

xyz 

Time 

(Hr) 

1 

000 

0 

41 

182 

18 

2 

010 

.01 

42 

202 

20 

3 

020 

.  02 

43 

222 

22 

4 

030 

.03 

44 

242 

24 

5 

040 

.04 

45 

262 

26 

6 

0  50 

.05 

46 

282 

28 

7 

060 

.06 

47 

302 

30 

8 

070 

.07 

48 

352 

35 

9 

080 

.08 

49 

402 

40 

10 

090 

.09 

50 

502 

50 

11 

100 

.  10 

51 

752 

75 

12 

150 

.  15 

52 

103 

100 

13 

200 

.  20 

53 

153 

150 

14 

250 

.  25 

54 

203 

200 

15 

300 

.30 

55 

253 

250 

16 

400 

.40 

56 

303 

300 

17 

500 

.  50 

57 

403 

400 

18 

600 

,60 

58 

503 

500 

19 

700 

.70 

59 

753 

750 

20 

7  50 

.75 

60 

104 

1000 

21 

800 

,80 

61 

154 

1,  500 

22 

900 

.90 

62 

204 

2,  000 

23 

101 

1.0 

63 

304 

3,000 

24 

151 

1.  5 

64 

404 

4,000 

25 

201 

2 

65 

504 

5,000 

26 

301 

3 

66 

754 

7,  500 

27 

401 

4 

67 

105 

10, 000 

28 

501 

5 

68 

255 

25,  000 

29 

601 

6 

69 

505 

50, 000 

30 

701 

7 

70 

755 

75,  000 

31 

801 

8 

71 

106 

100,000 

32 

901 

9 

72 

256 

250,000 

33 

102 

10 

73 

506 

500,000 

34 

112 

11 

74 

7  56 

7  50,000 

35 

122 

12 

75 

107 

106 

36 

132 

13 

76 

257 

2.  5x10^ 

37 

142 

14 

77 

507 

5.0x106 

38 

152 

15 

78 

757 

7.  5x10^ 

39 

162 

16 

79 

108 

lo”^ 

40 

172 

17 

80 

999 

largest 
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Table  II- 2  Continued. 


K 

81 

170 

0.  17 

82 

330 

.33 

83 

840 

.84 

84 

362 

36 

85 

482 

48 

86 

722 

72 

87 

962 

96 

88 

123 

120 

89 

143 

144 

90 

173 

168 

91 

243 

240 

92 

343 

336 

93 

503 

504 

94 

723 

720 

95 

104 

1,030 

K 

xyz 

Time 

(Hr) 

96 

144 

1,  440 

97 

224 

2,  160 

98 

444 

4,  380 

99 

764 

7,  570 

100 

884 

8,  760 

101 

135 

13, 140 

102 

185 

17,  520 

103 

205 

20, 280 

104 

355 

35,040 

105 

445 

43,800 

106 

885 

87,  600 

107 

886 

876,000 

108  through  128 

•  Spaie, 

set  000  at  present. 
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